Biased Percolation on Scale-free Networks 
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Biased (degree-dependent) percolation was recently shown to provide new strategies for turning 
robust networks fragile and vice versa. Here we present more detailed results for biased edge 
percolation on scale-free networks. We assume a network in which the probability for an edge 
between nodes i and j to be retained is proportional to (kikj)~ a with fc, and kj the degrees of 
the nodes. We discuss two methods of network reconstruction, sequential and simultaneous, and 
investigate their properties by analytical and numerical means. The system is examined away from 
the percolation transition, where the size of the giant cluster is obtained, and close to the transition, 
where nonuniversal critical exponents are extracted using the generating functions method. The 
theory is found to agree quite well with simulations. By introducing an extension of the Fortuin- 
Kasteleyn construction, we find that biased percolation is well described by the q — » 1 limit of the 
g-state Potts model with inhomogeneous couplings. 

PACS numbers: 



I. INTRODUCTION 

In recent years, much attention has been devoted to 
the study of real-life networks. Such networks may be 
modelled by points or nodes connected by edges. One 
feature is the scale-free topology, described by a proba- 
bility distribution P(k) for the number of edges k of a 
node, which falls off as a power law k~ 7 for large values 
of k. Most of the investigated cases turn out to have 
a topological exponent, or "degree exponent", 7, in the 
range 2 < 7 < 3.5. For 7 > 2 the mean degree (k) is 
finite, and for 7 > 3 also the variance (k 2 ) is finite. Like 
fully random Poisson-distributed networks (with a typ- 
ical scale), also scale- free networks are of "small- world" 
type. By now, many properties have been revealed and 
investigated thoroughly: these include degree-degree cor- 
relations, clustering and directedness of the edges in the 
network [H, 0, H, H • 

Another well-known property of scale-free networks 
is their resilience against random failure, a robustness 
caused by the presence of hubs (nodes with very high 
degree). On the other hand, these hubs may cause the 
network to be very vulnerable when a targeted attack is 
performed. In the limit of infinitely large networks, the 
network is said to be robust when even after removing 
an arbitrary fraction of the edges, there is still a nonzero 
probability that two randomly chosen nodes are part of 
a connected cluster. On the other hand, when remov- 
ing edges from a fragile network, a point will be reached 
when the giant cluster, the one with a size comparable 
to the network size, is destroyed; this very point is called 
the percolation threshold. The percolation transition is a 
genuine phase transition and is normally of second order 
so that critical exponents can be properly defined [ESQ- 

Since the first studies of percolation on scale-free net- 
works @, a lot of work has been done on node percola- 



tion 0, g |, 0, HH, O , bond percolationJTT|JTlJp, [H[ , 
percolation on multitype networks [H, Il4 Il5l| l6ll. clus- 
tered networks [l^ . Jl8l . correlated networks pL \& 1 1 7l [hH . 
directed networ ks M Jlrjl Il9j , degree-dependent edge per- 
colation [l3|, [3 l2(j] and degree-dependent node percola- 
tion [13 . 

The percolation transition has many connections to 
real systems. For example, it can be related to disease 
propagation models [H, [U [22| , [59[ . In this analogy, the 
infection of an individual is represented by the activation 
of a node of the (social) network. When a giant clus- 
ter of active nodes emerges, an epidemic is established. 
Disease propagation on such networks can be efficiently 
suppressed by selective vaccination, depending for exam- 
ple on the connectedness of each node. 

An alternative interpretation of a network with a cer- 
tain fraction of deactivated edges is in terms of a trans- 
port network in which the edges transmit data or deliver- 
ables between nodes with a certain transmission probabil- 
ity. This probability depends in general on the degrees of 
the connected nodes. For example, communication with 
the highly connected hubs on the internet is in general 
more efficient. It is, however, also possible that nodes 
with more edges are less robust. Indeed, in more so- 
cial terms, friendships involving people which have many 
acquaintances are more likely to end than friendships be- 
tween people with few connections. Or, as another ex- 
ample, traffic on a network induces high loads on highly 
connected nodes which in turn makes them more vulner- 
able to failure. Clearly, the resilience of an edge in a real 
network may depend strongly on the degrees of the nodes 
it connects. 

We study the properties of a network after biased or 
degree-dependent edge removal. More specifically, we 
consider networks in which the edge between nodes i and 
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j is retained with a probability proportional to its weight 



W; 



(kikj) 



(1) 



where ki and fcj denote the degrees of nodes i and j 
respectively, and a is the "bias exponent" . By tuning 
a, we can explore three qualitatively different regimes: 
random failure (a = 0), the attack of edges connected 
to hubs (a > 0) , and the depreciation of edges between 
the least connected nodes (a < 0). Henceforth, we call 
the regime a > "centrally biased" (CB). The converse 
regime, a < 0, is termed "peripherally biased" (PB). 

A degree dependence similar to that in Eq. ([T|) has 
already been considered in Refs. and [24] where Ising 
spin couplings Jy on scale-free networks were taken to 
be proportional to Wy . The motivation for introduc- 
ing degree-dependent couplings was the observation that 
for 7 < 3 the system is always "ordered" (critical tem- 
perature T c — oo) due to the dominance of the hubs. 
However, degree-dependent couplings make it possible to 
compensate high degree with weak interaction (assuming 
a > 0) so that the effect of the hubs can be neutralized. 
In doing so, it was discovered that a network with "inter- 
action exponent" a and degree exponent 7 has the same 
critical behavior as a network with interaction exponent 
zero (uniform couplings J) and degree exponent 



7 — a 
I- a 



(2) 



In this way it was possible to "trade interactions for 
topology" and study the rich mean-field critical be- 
haviour, with nonuniversal critical exponents depending 
on 7 ^, simply by varying a in a given network with 
fixed 7. The same exponent mapping will be recovered 
in this work in the following sense: at percolation the 
properties of a network with bias exponent a and degree 
exponent 7 are the same as those of a network with bias 
exponent zero and degree exponent 7, or degree exponent 
7, depending on conditions that will be specified. 

The significance and potential usefulness of biased de- 
preciation of a network is now becoming more clear. In- 
deed, it has been shown that networks with 7 > 3 are 
fragile under random failure, while networks with 7 < 3 
are robust under random removal of edges or nodes [1]. 
If it should turn out, and under certain conditions this 
is what we find, that the depreciated network behaves 
as one in which 7 is replaced by 7, it becomes possi- 
ble to control the robustness or fragility of a network 
systematically by tuning the bias exponent a. In other 
words, a network that is robust under random failure 
may turn out to be fragile under biased failure, and the 
other way round. Note that applying bias does not pre- 
suppose global knowledge about the network (location 
of the hubs, ...) but only requires local information on 
nodes and their degree. 

The exponent equality can be intuitively understood 
from the following heuristic argument, which is safe to 
use provided a > and k is sufficiently large. Using 



Eq. (fTJ) , one can anticipate that after depreciation of the 
network, a node with degree k will, on average, have a 
new degree k proportional to fc 1- ". Since all nodes re- 
main in place during the depreciation process, the origi- 
nal degree distribution P(k) changes into a new distribu- 
tion P(k) after depreciation, the relation between them 
being: 



P(k)dk = P{k)dk. 



(3) 



Using k oc k 1 a , one directly infers that indeed P(k) oc 

k 7 and the network after depreciation thus acquires de- 
gree exponent 7 and the corresponding percolation prop- 
erties. A more rigorous proof of this plausible expecta- 
tion is given in the Appendix. 

The paper is organized as follows: In Sees. Ill Al and lll B1 
we introduce random scale-free networks and present two 
distinct approaches by means of which a network can be 
reconstructed in a degree-dependent manner. Based on 
these schemes, we focus in Sec. Ill Cl on the degree distri- 
bution and the degree-degree correlations of the network 
after (partial) reconstruction. The percolation thresh- 
old is then extracted from these degree characteristics in 
Sec. IIIII The theory of generating functions for degree- 
dependent percolation on random networks will be exten- 
sively presented in Sec.[lV] In Sec. |V] the equivalence of 
our model with the Potts model is elaborated and using 
this equivalence and finite-size scaling theory, we arrive 
at the critical exponents for the percolation transition in 
Sec. I VII Finally, our results are extensively compared to 
simulational results in Sec. IVIII Our conclusions are pre- 
sented in Sect. I Villi A summary of part of the results 
presented here has been reported in Ref. [25|. 



II. DEGREE-DEPENDENT PERCOLATION ON 
RANDOM GRAPHS 

This section concerns (maximally) random scale-free 
networks. These are networks generated with the so- 
called configuration model, which assumes that the de- 
grees of the nodes in the network are distributed accord- 
ing to a probability P{k) which is taken to be the power 
law: 



P(k) = Ck~~<, 



(4) 



for values of k between the minimal and maximal degrees 
m and K, respectively. C is the normalization constant. 
In order to ensure a finite mean degree we take 7 > 2. 
The graph is then completed by connecting the stubs 
emanating from all nodes. The probability P n {k) that a 
randomly chosen edge leads to a node of degree k must 
therefore be: 



Pn(k) 



kP(k) 



(5) 



where (•) denotes the average over the nodes, obtained 
using probability distribution P(k). The probability dis- 
tribution P n is also called the nearest-neighbor degree 
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distribution. Random networks are constructed by con- 
necting the earlier mentioned stubs randomly. We do not 
allow self-connections nor multiple connections between 
nodes and use the method proposed in Ref. [44] to avoid 
degree-degree correlations in the network. To quantify 
degree-degree correlations, let us introduce the probabil- 
ity P(k, q) that nodes of degree k and q are connected. 
If no correlations are present, P(k,q) reduces to: 



P(k,q)=P n (k)P n (q) = 



kqP(k)P(q) 



(6) 



Below and close to the critical point, large and random 
networks can locally be treated as trees and loops are 
sparse so that their effect can, to a good approximation, 
be ignored. The local tree-like structure will be used in 
Sec. IIVI when the generating functions method is intro- 
duced. 

We continue with presenting two distinct depreciation 
methods to study degree-dependent percolation. The 
statistical edge properties are now being considered, 
which will allow us in Sec. Ill CI to obtain the statistics 
of nodal properties. 



A. Sequential Approach 

Our first method, which we call the sequential ap- 
proach, starts from a random network with all N node 
degrees distributed according to the degree distribution 
P(k). Initially all edges are removed and we aim at 
reintroducing a fraction / of the total number of edges 
N e = (k)N/2. This is achieved by activating one edge 
in each time step t. Consequently, the probability that 
the edge between nodes i and j is activated is Wij/Z t 
where Z t is the sum of weights Wij of all non-activated 
edges after t— 1 steps. Thus, the probability Pij(f) that 
an edge between nodes i and j is again present after the 
reinclusion of a fraction / of the edges, is [26| : 



Pij(f) 



n 

t=i 



W; 



(J) 



For sufficiently large networks Wij/Z t is typically small 
compared to one and Eq. Q is well approximated by: 



putt) 



-DfWij 



(8) 



with the positive parameter Df = Ylt=i ^t 1 ■ It can be 
argued that for a sufficiently narrow distribution of the 
weights foOj : 



Z t = (w) e {N e -t + l), 



(9) 



where (-) e denotes the average over all edges. The fol- 
lowing property is readily derived, 



7-2 



V{w) ^^kT l-2 + a 



(10) 



Using Eq. ([9]), Df can be determined, such that for 
large N e , 



and thus W£ 



D f = -]n[l-f]/(w) e 



Pij (f) = !-[!- /]«**/<»>• 



(11) 



(12) 



It is instructive to consider a few asymptotic regimes of 
Eq. (fT2")l . First, in the case a = 0, one recovers the 
expression for degree-independent percolation pij = f as 
expected. Secondly, for arbitrary a, we can distinguish 
the dilute limit and the dense limit in terms of /, and 
find: 



Pij ~ fWij/{w) e when / — > 0, 
Pij ~ 1 when / — > 1. 

Thirdly, when w lj /(w) e < [- In (1 - 

Pij ~ -ln(l- f)w lj /(w) e . 



(13a) 
(13b) 



(14) 



We proceed by defining the marginal distribution pk as 
the mean probability that an edge connected to a node 
with degree k is present in the network after reconstruc- 
tion. Thus 



K 



Pk 



E p «(«)a 



kq > 



(15) 



q—7n 



where pk q = 1 — e ~~ D f Wkc '. A good analytic approxi- 
mation to pk can be obtained by substituting p^ q into 
Eq. (fl"5)) . considering an integral instead of a sum, tak- 
ing the macroscopic limit K — > oo, and expanding the 
exponential, 



Pk = 1 - 



n=0 



(-)"(£>/ wC a k- a ) n 
(1 +na/{j- 2))n\ 



(16) 



Alternatively, this result follows straightforwardly from 
the fact that the marginal distribution involves the in- 
complete Gamma function [46]. The usefulness of this 
explicit form can best be appreciated by first considering 
the range < a <C 7 — 2, for which we obtain the simple 
analytic result 



Pk 



1 — exp 



7 



D f mT a k- 



= l-cxp{-D f ^)~ e k- a } 



(17) 



Although this result is strictly only valid for the specified 
range of a specified above, numerical inspection shows 
that it is a rather good approximation to the integral 
representation of the sum Eq. (I15p for a wider range of 
a, including negative values. In fact, the result is useful 
in the entire interval of our interest a G [2 — 7, 1]. Using 
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the previously obtained approximation to D/, Eq. (fTTj) . 
it can be further simplified to 



Pk 



i - [i - f) k a 'V<vT' 



(18) 



Based on the asymptotic regimes of ptj it is also possi- 
ble to extract the behavior of p k - When / — > 0, we may 
use Eq. (fl3a|) : 



Pk ~ k- a f/^{w) ( 



(19) 



On the other hand, when a > and fc S> /c x , where the 
cross-over value for k is given by k x = dV" 1 /m, we get: 



Pk ~k- a D f ^/tw) t 



In(l-/)/V(«5)e. (20) 



B. Simultaneous Approach 

As an alternative to the sequential approach we in- 
troduce now the simultaneous approach. Again we start 
from a fully depreciated uncorrelated network with de- 
gree distribution P(k). We then visit each edge (between 
nodes i and j) once and activate this edge with proba- 
bility 



Pij = fWij/(w) e . 



In contrast to the sequential approach, p^ is now history- 
independent. Note also that (pij) e = / as it must be. For 
the marginal distribution p k one finds: 



Pk = k- a f/^w)~e, 



(22) 



which satisfies fpkq = PkPq and is the same as in the / — > 
limit of the sequential approach (Eq. (|T9]) ) . However, 
for each value of k and q, the probability p kq must be less 
than or equal to one. This means that Eq. (|2"Tj) is only 
well-defined for values of / for which 



/ < fu 



(A- 1 



(k) 



x (min(m Q , K a )f 



(23) 



It can be calculated that, in the macroscopic limit, the 
rhs of Eq. (|23|) vanishes when a < and therefore the 
simultaneous approach is only meaningful for positive a 
and provided 



/< 



7- 



7 



(24) 



To reach fractions above this limit in the simulations, 
we iterate the simultaneous approach. The first iteration 
involves the usual simultaneous approach with / = f u ; 
the second iteration is initialized by considering a new 
network consisting of all edges that have not been rein- 
troduced during the first sweep. For that network one 
calculates the probabilities Wij/Z 2 and a new value of 
/„, which is the minimum of the set {{w) e /wij} where 



also the average is only over edges of the new network. 
One then applies the simultaneous approach until the 
new f u is reached, after which a third iteration can be 
initialized if necessary. Such iterations, however, intro- 
duce correlations and history dependence. Note that the 
sequential approach can be seen as an extreme case of an 
iterated simultaneous approach in which only one edge is 
reconstructed in each iteration. 



C. Degree Distribution and Correlations of the 
Reconstructed Network 

We now seek to obtain the degree distribution and 
characterize degree-degree correlations for the network 
after reconstruction. The following is valid for both the 
simultaneous and sequential approaches. Henceforth we 
adopt the convention that an overbar indicates quantities 
in the diluted, or depreciated, network. 

For the node degree distribution P(k) and the degree- 
degree correlations embodied in P(k,q) of the depreci- 
ated network we can write 

K 



(21) 75.T 



P(k)=^ L P(A k AB k ^), (25a) 

k=k 
K K 

P(k,q) = YJ1_ P{ * C ^ A A B ^k A D )/f- (25b) 

k=k 9=9 

Here we introduced the notation, for events A-D, 

• A k : a randomly chosen node of the original net- 
work has degree k. 

• Bk^k : tlie degree of a node goes from k in the 
original to k in the depreciated network. 

• C q k'. the nodes connected by a randomly chosen 
edge of the original network have degrees q and k. 

• D: the chosen edge has not been removed from the 
original network. 

For the node degree distribution, one readily finds 



K 



p{k) = Y J m 



ft-* 



(26) 



For large values of k, i.e., k 3> k x , and a > 0, the proba- 
bility of retaining a node of degree k falls off as p k oc k~ a \ 
this is valid using the sequential approach (see Eq. (1201) ). 
as well as the simultaneous one (see Eq. (|2"2"10 . Substi- 
tuting this into Eq. (|26[) and approximating the binomial 
distribution in Eq. (|26p by a normal distribution, one 
arrives at: 



P(k) oc k for k 



(27) 



where 7 is defined in Eq. ([3]). This result, which is proven 
in the Appendix, confirms the validity of the expectation 



5 



raised in the Introduction. Note that in case a = 0, 
P(k) oc fc 7 as it must be. We introduce now averaging 
over nodes of the reconstructed network: 



K 



< ■ >= 



E p ( fc ) 



(28) 



fe=0 



For further purposes, we calculate now the first and sec- 
ond moment of P(k) in terms of the moments of P(k): 



«fc» = (fcp fe )=/(fc), (29a) 

<fc 2 > = (kpk(k Pk - pk + l)). (29b) 

The degree-degree correlations are embodied in 
Eq. (|25b[) . This function can be further worked out to 
yield: 

K K 

P(k,q) =^p n (9)Pn(%fe, 

i j-zDdTHi-Pk)"-*//. (so) 



k=k 1=1 
X 



This can be reduced to: 



K K 



P(k,q) -fcg^E 



k=k 1=1 



(f(k)) 2 



PkPq 

(31) 



Eq. (|3"Tj) expresses the degree-degree correlations of a net- 
work after degree-dependent depreciation of a fully un- 
corrected network. The question that can now be raised 
is when the depreciated network is also free of correla- 
tions, or, when is P(k,q) = P n (k)P n (q)7 It is readily 
checked that this is true provided 



fPkq = PkPq- 



(32) 



As Eq. (f3"2")l is valid for the simultaneous approach (see 
Eqs. (|2"2"|) and (pH]) ). no correlations appear in the re- 
constructed network (after a single iteration). For the 
sequential approach, on the other hand, Eq. (|32p is gen- 
erally not satisfied and the reconstructed network will be 
correlated. 

The following limiting case of the sequential approach 
is interesting: take a positive and consider an edge be- 
tween two nodes of large degrees k and q such that 
k >• fc x and g> k x . We may then substitute Eqs. (fT3|) 
and (|20p into Eq. (f3Tj) . One soon arrives at the result: 



P(k,q) 



f 



Pn(k)P n (q) In (I-/)' 



(33) 



Since the rhs is smaller than one, this demonstrates that 
the sequential approach causes disassortative mixing in 



the depreciated network when a > 0. In other words, 
nodes with large degrees tend to be connected to nodes 
with small degrees and vice versa. Using simulations, we 
will present evidence in Sect. IVIll that such correlations 
are introduced. 

Finally, note also that Eq. ([3Tj) reduces to the correct 
nearest-neighbor degree distribution upon summing over 

q- 



Pn{k) 



kP(k) 



(34) 



III. 



PERCOLATION THRESHOLD FOR 
CENTRAL BIAS 



Here we focus solely on centrally biased (CB) deprecia- 
tion (a > 0) using the simultaneous approach. One may 
wonder what happens if centrally biased (CB) deprecia- 
tion is applied to a robust network with 7 < 3 such that 
the edges between and emanating from hubs are prefer- 
entially removed. Since, in that case, 7 > 7, one may 
speculate that a robust network may turn fragile and 
that the threshold for this to occur is 7 = 3 instead of 
the threshold 7 = 3 valid for degree- independent percola- 
tion. We will address this question further and conclude 
that it is indeed so. 

On the other hand, if we start from a fragile network 
(7 > 3) and apply CB, it is logical that the net remains 
fragile. Upon removing edges linked to hubs with a larger 
probability, we are more likely to destroy the coherence 
of the network. The question can then still be posed 
how much the percolation threshold of the reconstructed 
network is shifted. 

Our first task now is to calculate the critical fraction 
at which the network becomes disconnected. According 
to Molloy and Reed [43j , the critical fraction of a random 
network can be found by looking at the average nearest- 
neighbor distribution. If, upon following a random edge, 
the attained node has more than two neighbors, the net- 
work is said to be percolating, that is, a giant cluster 
will be present in the network. Note that this criterion 
is exact for the simultaneous but not for the sequential 
approach, due to the appearance of degree-degree corre- 
lations. In the reconstructed network, the Molloy-Reed 
criterion reads: 



< k{k - 1) > 



(35) 



or e 



« k > 
;quivalently, using Eq. (|29|) : 

2(kp k ) = (kp k (k Pk - p k + l)). (36) 

Using Eq. (|22"|) for the simultaneous approach, we find the 
following expression for the critical fraction f c at perco- 
lation 125]: 



fc = 



(k 1 



(k)((k 



2~2a\ 



(fc 1 - 2 "))' 



(37) 
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For unbiased depreciation of the network (a — 0) this last 
expression reduces to the well-known formula of random 
percolation on random networks Q. 

Eq. (|37[) allows us now to find out whether the network 
is robust, or in other words, whether f c — ► 0. This van- 
ishing occurs, in the macroscopic limit, when the term 
(k 2 ~ 2a ) diverges and therefore: 



> 3: the network is fragile, 
< 3: the network is robust. 



(38) 



The scaling relation of f c as a function of the network 
size can be found when 7 < 3; using N oc K 7 ^ 1 , one 
finds: 



(39) 



Whether or not a network is robust for the degree- 
dependent attack is thus not solely a property of the net- 
work. Also the exponent a plays a crucial role in the 
arguments and its effect can be absorbed by using the 
exponent 7 instead of the exponent before dilution, 7. In 
Sect. IVII we will take a closer look at the regime around 
the percolation threshold and we will find that the same 
mapping from 7 to 7 is valid. 



IV. GENERATING FUNCTIONS APPROACH 

We now introduce the generating functions approach 
for degree-dependent percolation. By this method, cer- 
tain properties of the finite clusters in the network are 
easily obtained; this in turn allows to draw conclusions 
about the giant cluster. The method is exact if loops in 
the network can be ignored; since in the macroscopic 
limit, the average loop sizes in the finite clusters di- 
verge [57j . the method turns out to be exact. This will 
be apparent in Sect. IVTil when comparing the analytical 
results with simulations. 



We also define the generating function for the distribu- 
tion of residual edges of a node reached upon following a 
random edge: 



K 



-h(k-l) 



(41) 



The exponent of e~ h is k — 1 because the edge which is 
used to reach the node is not counted. 

There is a threefold advantage in working with the 
generating functions instead of working with the degree 
distribution itself. First, moments can be obtained easily 
from the generating functions. For instance, the average 
degree is given by 



K 



(k) = kP{k) = -G' (0), 



(42) 



where G' denotes the derivative with respect to h. 
Higher moments can be obtained with higher-order 
derivatives. 

Secondly, we can benefit from the so-called powers 
property of generating functions: if the distribution of a 
property k of an object is generated by a function G(h), 
then the generating function of the sum of n independent 
realizations of k is G(h) n . For instance, if we randomly 
choose n nodes in our network, the distribution of the 
sum of the degrees is generated by Go(h) n . 

Thirdly, the use of generating functions will allow us 
in Sect. [V] to highlight the equivalence with the q — > 1 
limit of the g-state Potts model where the parameter h 
will play the role of the magnetic field. 



B. Self-consistent equations 

We can now define the equivalents of Go(h) and G\(h) 
for the network after dilution as Fq[K) and F\{h): 



A. Introduction 

Generating functions are used in a wide branch of 
mathematical problems concerning series (35j |. A gen- 
erating function of a series is the power series which has 
as coefficients the elements of the series. Applied to the 
context of percolation problems, this series is taken to 
be that of the discrete probability distributions charac- 
terizing the network under consideration [IH, |47| . We 
explain first the general formalism while closely follow- 
ing the approach of Newman [l3j |. which we adapt for 
degree-dependent edge percolation [6l| . 

The most fundamental generating function is the one 
that generates the degree distribution of the network 



-hk 



K 



G (h) = P(k)e 



-hk 



(40) 



F (ft)=x;p(fc)t 

fc=0 
K 

F 1 (h) = Y / Pn(k)e- h(k - 1) . 



(43a) 
(43b) 



fe=0 



Note that the minimal degree in the network after dilu- 
tion is zero instead of m. Eq. (|4"3"]) can be worked out 
further using Eqs. ([26]) and ([331): 



A" 



F (h)= J2 p ( k )( l ~Pk + ^ h Pk) k , 



(44a) 



K 



F 1 (h) = ^^(1 - Pk + e- h p k ) k -\ I 441, i 



k=m 



f 



The most interesting quantity for us is the size distri- 
bution of the finite clusters, the generating function of 
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which can be readily derived using Fq and F\. Let Hq 
denote the generating function for the probability that 
a randomly chosen node belongs to a connected cluster 
of a given (finite) size. Furthermore, let Hi be the gen- 
erating function for the probability that upon following 
a randomly chosen edge to one end, a cluster of a given 
(finite) size is reached. If the network can be treated 
as a tree, these generating functions satisfy the following 
self-consistency equations [62j : 



C. Full Derivation of Self-consistent Equations 



H x {h) = e- h F 1 [H l {h)], 
H (h) = e^FolH^h)]. 



(45a) 
(45b) 



Here the function Fo } i[Hi(h)] denotes the function Fo,i 
in which e~ h is replaced by H\(h). The proof of these 
relations relies on the aforementioned powers property 
and is expounded in Ref. [l3|. The percolation threshold 
can now be derived with the aid of these functions. 

Several macroscopic quantities can be easily identified 
in the depreciated network Q. For example, we define 
Too as the probability that a node belongs to the giant 
cluster, Coo as the edge probability for being in the giant 
cluster and S as the average cluster size of finite clusters: 



We prove now that the self-consistent Eqs. ([45]) are 
only valid in case no correlations are introduced in the 
reconstructed network, or, when fpkq = PkPq, as is valid 
for the simultaneous method only. Here we will give a 
precise derivation of the self-consistent Eqs. ([45]) . thereby 
taking into account the degree-dependence of the func- 
tions Hi. 

Let us first look at the generating function Hf~^ (h) 
for the probability that an edge, which connects nodes of 
degree q and k, branches out in a cluster of a given edge 
number along the node of degree k. It is readily derived 
that Hf^ k satisfies the equation: 



K 



Hr\h) = e -Ml+£ P n {k)p k% [H^ k {h) 1] 



k—m 



k-1 



(51) 



H*{h) 



We proceed by defining 
Y,k P ni. k )PqkHl~^ k {h)/ p q , such that 
set of self-consistent equations for each value of q: 



we arrive at a 



Voc = l-ffo(O), 
£ oo = l-(i/ 1 (0)) 2 , 
S = -fl£(0). 



(46a) 
(46b) 
(46c) 



K 



Moreover, the degree distribution of nodes in the giant 
cluster varies as 



P gc (k) «(!-(#! (0)) fe )P(fc) 



(47) 



from which it follows that the degree distribution of 
the finite clusters varies as Pf c (k) oc (Hi(0)) k P(k). In 
case there are both finite clusters and a giant cluster, 
the degree distributions have the asymptotic behavior 
(k — > oo): 



P 9C (k) ■ 
Pfc(k) 



-fc/A 



(48a) 
(48b) 



with A = — ln(i?i(0)). In other words, in the presence of 
a giant cluster, only the degree distribution of the giant 
cluster falls off with a power law with exponent 7. The 
average cluster size (I46cp in the diluted network, on the 
other hand, can be further worked out by differentiating 
Eqs. ([4"5]) with respect to x: 



5=1 



f(k) 



1 + F>(0) 

Hence the average cluster size diverges when 



1 



-F[(0). 



(49) 



(50) 



This is yet another way of writing the Molloy-Reed cri- 
terium Eq. 



Hf(h)=e- h Y / 



Pn(k)p q k 
Pq 



l + Pk [H k (h)-l]) 



k-l 



(52) 



Note that, as derived in Sec. Ill CI no correlations are 
induced during depreciation when f p kq = PkPq- In that 
case, this equation reduces to Eq. (j45a|. After solving 
Eq. (fS"2"|) with respect to H k for all values of k, we can 
also calculate: 



K 



H (h) = e ~' lp (fc) [1 + Pk[H k (h) - 1] 



(53) 



Again, this expression reduces to Eq. (|45b() in case of 
correlation-free depreciation when H\ is independent of 
k. 

Below the percolation transition, a trivial solution ex- 
ists: H k (0) = 1. This solution, however, turns unstable 
at the percolation threshold. The threshold value may be 
derived by linearization of H\ (0) around its equilibrium 
value: H\ (0) = 1 — Ek with e k <C 1. The percolation 
criterion is then: 



K 



yy P n (k)pqkPk ^ 



Pq 



l)e k - 



(54) 



Again, in absence of correlations in the network, that 
is, when the criterion f p q k = PkPq is satisfied, this re- 
duces to the earlier encountered Molloy-Reed criterion 
of Eq. ([36]) . Eq. ([54]) is the criterion for the percolation 
threshold for correlated systems, such as the one created 
using the sequential method. However, solving Eq. (|54[> 
to obtain f c constitutes a rather difficult task. 
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D. Original and Depreciated Network 

We show now that another approach exists by which 
one easily derives the self-consistent equations charac- 
terizing the network. This method is closer to the one 
followed in other works. 

Let us call the probability that an edge in the 

network does not lead to a vertex connected via the re- 
maining edges to the giant component (infinite cluster) 
and pij the probability that the edge between nodes i 
and j is active. Then, following the edge along node j, 
one finds: 



R^ = 1 - Pij + Pij f[ B?- 



(55) 



z — l...kj 



This type of equation was already obtained for the q — > 
1 limit of the g-state Potts model [Hj]. Using the tree 
approximation, we can rewrite everything as a function 
of the degrees of the nodes: 



R" 



-^k 



1 — Pqk + Pqk 



K 



P n(k)R k - 



k-l 



(56) 



This self-consistent set of equations is equivalent to the 
ones that we obtained in Eq. ([52)1 . Indeed, after the 
transformation R q ^ k — 1 = p q k(R q ^ k — 1), we find: 



R q - 



K 



E Mk) Pkk {R k - 



k—rn 



k-l 



(57) 



which is exactly the same as Eq. (f5T|) for h = 0. 

The purpose of this derivation is to show that our self- 
consistent equation (|5ip is in agreement with Eq. (|56[) . 
For degree-independent R and p, an equation similar to 
Eq. ([56]) appears frequently in the literature. The differ- 
ence between Eq. (|51l) and Eq. ([57]) is that Hi is normal- 
ized with respect to the depreciated network, whereas R 
is normalized with respect to the original network. 



V. EQUIVALENCE WITH THE POTTS MODEL 

There exists an equivalence between edge percolation 
and the q — + 1 limit of the g-state Potts model. This 
connection was first worked out by Fortuin and Kaste- 
leyn in Ref. [H. Although initially used for lattice mod- 
els, the connection was very general and is valid for any 
network [2(| [38J. Moreover, their proof can easily be 
generalized to incorporate edge-dependent coupling con- 
stants and edge-dependent removal into the Potts model 
and the percolation model, respectively. We explain here 
in more detail this equivalence and reformulate our per- 
colation problem as a spin-like problem which will allow 
us to derive critical exponents and compare them with 
the ones obtained for the Potts model. We will also find 



support for our simple scaling relation using exponent 
7 (Eq. ©), as was already encountered in studies con- 
cerning degree- dependent Ising interactions on scale-free 
networks Hi, |24|. Note that the Ising model [U 53 
and the Potts model [HI, 42] , together with their critical 
properties, were already studied on scale- free networks. 

The Potts model can be seen as a generalization of 
the Ising model in which each site i has a spin ui. In 
the Potts model, these spins can take q distinct values 
0, . . . , q — 1 and the Potts Hamiltonian is: 



(58) 



Here < ij > indicates nearest- neighbor sites i and j, Jij is 
the coupling constant and 5 the Kronecker delta function. 
Note that the Ising model corresponds to the q = 2 Potts 
model. The Fortuin-Kasteleyn theorem states now that 
the free energy of the q — ► 1 limit of the g-state Potts 
model is the same as the "free energy" of the percolating 
network where the latter is the generating function of the 
cluster size distribution function n,: 



^ — hs 



(59) 



Here the average is performed over all networks in which 
the probability to retain the edge between nodes i and 
j is pij. The parameter pij in the percolation problem 
corresponds in the following way to parameters of the 
Potts model M [H: 



Pij 



1 



-Jij/k B T 



From this relation, we can immediately identify the prob- 
ability Too for a node to be in the infinite cluster and the 
average cluster size S, earlier introduced in Eq. (jlu) . As 
we are interested in the behavior near criticality, we in- 
troduce 



e = f-f c 



and obtain 



Vooie) = 1 + 



dT_ 
~dh 



h=0 



S(e) 



d 2 F 



dh 2 



(60) 

(61a) 
(61b) 



h=0 



Note also that J-(e, 0) gives the total number of finite 
clusters. 



VI. SCALING THEORY AND CRITICAL 
EXPONENTS 

In the following section, we introduce finite-size scal- 
ing in order to find critical exponents near the percolation 
transition. In order to solve the scaling relation, we use 
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a Landau-like theory which we derive from the exact re- 
lations (l4"5l) . We follow closely the approach presented in 
Refs. l53l l54l andU^ for finite-size scaling in systems with 
dimensions above the upper critical dimension. However, 
we will find that the forms of the Landau-like theories of 
Refs. IH, and[3(3 were too limited for studying the per- 
colation transition in case the distribution function has 
a very fat tail, that is when 2 < 7 < 3. 

According to finite-size scaling, the free energy J 7 of a 
large but finite network with N nodes close to criticality 
can be written in the general form [48| : 



T 



(e, h) = AT- 1 ! (eN l ' v % hN x ' Vh } ( 



(62) 



where F is a well-behaved function. The variable eN 1 ^' 
originates from the existence of a "correlation number" 
N% (instead of a correlation length) which scales as N% oc 
e~ u " such that the first variable of F can be rewritten 
as (AT/AT^) 1 /^ [5^ [jyj. it is then obvious that close to 
criticality, the (singular part of the) free energy scales as: 



T(e,0) oc e u % 
F(0,h) oc h Uh . 



(63a) 
(63b) 



As a second scaling ansatz, we assume that, in the macro- 
scopic limit, the scaling of the cluster size distribution n s 
can be written as [H| : 

n s (e) = s- T G(es ff ), (64) 

where again G is a well-behaved function. The form of n s 
one usually has in mind is n s oc s~ T e~ s / s |13j . which is 
essentially a damped power-law with cutoff s* , valid for 
large cluster sizes. At criticality very large clusters arise, 
caused by the diverging cutoff s* according to s* oc e -1 /' 7 
such that n s (0) ~ s~ r . 

Using the analogy with the Potts model, we define the 
usual critical exponents a, j p , (3 and 5 for the percolation 
problem as: 



^(e,0) 
^00 (e) ' 
5(e)- 

dT_ 

~dh 



+ 1 



(65a) 
(65b) 
(65c) 

h 1 ' 8 . (65d) 
and (|6"4")l , standard 



Using the scaling forms of Eqs. 
techniques provide us with exponent relations by which 
all critical exponents can be related to Vh and v t . One 
arrives at [48j 



13 = 




(66a) 


7p = 


v^h 1 - !). 


(66b) 


a = 


2-2/i-7p, 


(66c) 


a = 


(P + lr>)-\ 


(66d) 


T = 


2 + P((3 + lp r\ 


(66e) 


5 = 


(/3 + 7 P )/A 


(66f) 


V S = 


P + 1p- 


(66g) 



The last exponent Vf is related to the usual fractal di- 
mension df in the same way that v e is related to the usual 
dimension d and quantifies how the cluster size s scales 
with e close to criticality, i.e., s oc e~ Vf . Note that the 
cutoff size s* scales like s, since vj = 1/a. 

The problem we are left with now is to find the scaling 
exponents Vh and v t for percolation on scale-free net- 
works. This can be done in an exact way since we know 
the equation of state from Eq. (|45[) as a function of the 
order parameter: 



i/j(e,h) = l-H 1 {e,h). 



(67) 



As we are merely interested in the behavior near the tran- 
sition where e < 1, 4 < 1 and ip -C 1, we can expand 
Eq. (|45a|) . For the case 7 > 3, we find the form Q: 



h = -cieip + c 2 ^ 2 + ... + Csijj 1 2 + 



(68) 



in which all Ci as well as the coefficient of the singular 
term, c s , are positive constants. This equation also fol- 
lows from minimization of the free energy: 

T{e, h) oc -hip - ciap 2 + c 2 ip 3 + ... + c^ 7-1 + . . . , 

(69) 

with respect to the order parameter ip. 

We distinguish two cases now. First, when 4 < 7, 
we know that f c is finite and the relevant part of the 
equation of state for ip becomes: 



h = —cieip + C2V' 2 



(70) 



Solving for ip, and substitution into Eq. (|69[) one then 
simply finds that the free energy scales as: 



F(e,0) oc e 3 , 
F{0,h) oc h 3/2 . 



(71a) 
(71b) 



In other words, when 4 < 7, we find that v t = 3 and 
Vh = 3/2. From these two exponents, and using Eq. (|66|) . 
we list all other exponents in the last column of Table 1. 
Note that, as expected, these exponents a gree with the 
usual mean-field results for percolation |47l~l48l]. 

Secondly, when 3 < 7 < 4, the relevant part of Eq. (|68|) 
reduces to 



h = -cieip + c s tp J 2 , 

from which follows that 

.F(e,0)~eH 

F(0,h) - h%=*. 

Therefore, we come to: 

7-1 7-I 
Ve = - — - and v h = - — -. 

7—3 7^2 



(72) 

(73a) 
(73b) 

(74) 



and we obtain all other exponents as given in the second 
column of Table 1 . 



10 



TABLE I: Critical Exponents 





2<7<3 


3<7<4 


7 > 4 


13 


1 


1 


1 


3-7 


7-3 


r 


27-3 
7-2 


27-3 
7-2 


5/2 


a 


3-7 
7-2 


7-3 
7-2 


1/2 


a 


37-7 
3-7 


5 — 7 
7-3 


-1 


7p 


-1 


1 


1 


5 


7-2 


7-2 


2 




7-2 
3-7 


7-2 
7-3 


2 



Lastly, in case 2 < 7 < 3, the critical fraction / c = 0. 
Again, we expand Eq. (|45ap using the small parameters 
e <C 1, \h\ <C 1 and i/iCl The equation of state for ip 
and the associated free energy (f6"9"| become: 



h = -c s (e^- 2 + V + ••• 
where c, > 0. It follows that 



^/2) 



^(e,0) 



(75a) 
(75b) 

(76a) 
(76b) 



In the last expression, the limit e — > is only taken in 
the free energy and h is taken small and negative such 
that ijj is still positive. Therefore, 



7-I 7 - 1 

— and Vh = — 



3-7 



7-2 



(77) 



The other exponents are listed in the first column of 
Table 1. It must be noted here that in practice, the 
exponents of Eq. (|77|) may be impossible to find with 
the configurational model in case we start from a ro- 
bust network (2 < 7 < 3). This stems from the fact 
that a structural cutoff for the maximally allowed degree 
K must be introduced to obtain an uncorrelated net- 
work. Such cutoff can be of the form K ~ N 1 ^ with 
u> £ [2,oo[. However, it is well-known that the cutoff af- 
fects the critical exponents (58[. Indeed, performing the 
averages in Eq. (|3~7| with use of the cutoff, one readily 
obtains f c cx jv(i-«)(r-3)/w ( cf Eq an d therefore: 

(78) 



(l-o)(7-3) 



In case we start from a fragile network (7 > 3), u> equals 
7 — 1 and the z/ e of Eq. ([77)1 is retrieved. Therefore, 



degree-dependent edge removal may be used as a tool 
to observe critical exponents in the delicate regime 2 < 
7 < 3, without the use of a structural cut-off 58J]. In our 
simulations as presented in Sect. IVII1 we take w = 2 in 
case 7 < 3. 

The exponents obtained for the case 7 > 3 reduce to 
the ones for node percolation in the limit a — when 
7 = 7 [5j. However, in the limit a — the exponents 
in the first column do not coincide with those given in 
Refs.il and!. 

Close to the percolation transition, it is possible, with 
the use of Eqs. (|46|) to calculate the average degree in the 
giant cluster. We can identify this as 



/iV e (l-(^) 
/-»/. N{1 - Ho) 



(f c + e)(k)2N 
N(k)(f c + e) 



(79) 



This however, must not be confused with the Molloy- 
Reed Criterion, which states that the average degree of 
a neighboring site in the entire network has degree two. 

In sum, we have now calculated the most important 
critical exponents for a percolation process. Extra sup- 
port for our critical exponents comes from scaling rela- 
tions and the connection with the Potts model. 

There is one feature which appears in all the calcu- 
lated exponents: the only dependence on a arises through 
the exponent 7. Random percolation on a network with 
degree exponent 7 gives the same critical exponents as 
percolation with bias exponent a on a network with 
degree exponent 7. This equivalence was found be- 
fore for de gree -dependent interactions on scale-free net- 
works [23l |24j. It is not surprising that the same be- 
havior appears both for edge percolation and for degree- 
dependent interactions, since both can be linked with the 
Fortuin-Kasteleyn construction. 



VII. COMPARISON WITH NUMERICAL 
RESULTS 



In this section we test the previously derived analyti- 
cal results using simulations. The networks are generated 
using the uncorrelated configurational model which was 
introduced in Sect. [TTJ Each simulation involves three 
free parameters: the degree exponent 7 of the network, 
the minimal node degree m and the number of vertices 
N. Unless mentioned otherwise, we set m = 1. In the 
configuration model, the degrees of the nodes are deter- 
mined initially from the discrete degree distribution [65[ 
and then the connections are assigned at random. To 
obtain an uncorrelated network, the maximal degree is 
set to V^/V when 2 < 7 < 3 44 1 . In some simulations (see 
Figs. Q] and [2|) no degree cutoff was imposed. If 7 > 3, 
the maximal degree is simply N — 1. Both the sequential 
and the simultaneous approach are implemented. 
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FIG. 1: The critical fraction at percolation f c as a function 
of the network size N. To compute this critical fraction we 
average over 10 4 network realizations for each set of param- 
eters and apply for each network the percolation process 100 
times. No cutoff was introduced for the maximal node degree 
and network reconstruction was done with the sequential ap- 
proach. 
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FIG. 2: The probability p^ to retain a node after deprecia- 
tion (see Eq. (jT5j) ) as a function of its original degree k using 
the sequential approach. The main panel shows results for 
networks with 7 = 4 submitted to PB with a = — 1. The 
dots indicate the simulation results for ten network realiza- 
tions and ten percolation routines. The continuous lines are 
the best fit to the data of Eq. (|17|) . The inset shows the same 
but for a network with 7 = 2.5 subjected to CB with a = 0.5. 
In both cases, m = 2 and N = 10 s . The arrows indicate the 
crossover value k x ■ No cutoff was introduced for the maximal 
node degree. 



A. Sequential Approach 



In this first subsection, we discuss simulation results 
concerning the sequential approach. Most of these results 
can also be found in Ref. [251 . 



1. Scaling of the Critical Point 

First, we investigate the finite-size behavior of the crit- 
ical fraction f c of nodes. The results for the sequential 
approach are shown in Fig.[T] For networks with 7 = 2.5 
submitted to CB with an effective value 7 = 4 (contin- 
uous black line), we observe that the critical fraction f c 
converges to a finite value as N grows, confirming the 
conjecture that a robust network may turn fragile under 
CB. In the opposite network with 7 = 4 submit- 

ted to PB with an effective 7 = 2.5 (dashed red line), has 
a critical fraction that decays with the vertex number iV 
as a power-law, f c ~ N~ x / Vc . The best fit to the data 
in this case results in l/v e = 0.35 ± 0.02, consistent with 
the value 1/3 expected from Eq. l[3"5j) . This result shows 
that a fragile network under PB will behave in the same 
fashion as a robust network with a deg ree distribution 
controlled by 7 under random failure [66| . Note that this 
simulation result confirms Eq. ([55]) . although this equa- 
tion was derived for the simultaneous approach. Indeed, 
to deduce Eq. correlations in the diluted network 

are neglected. These simulation results indicate that this 
is an acceptable approximation. 



2. Properties of the Diluted Network 

During our theoretical discussion of the sequential ap- 
proach, certain properties of the diluted network became 
apparent, such as a cross-over behavior as a function of 
the degree k for the mean edge preservation probabil- 
ity pk- The probability pk can easily be inferred from 
our simulations by calculating the ratio of the new node 
degree k to the old node degree k for each node and av- 
eraging over all nodes with the same degree. The result 
is shown in Fig. [2] in which the expected crossover behav- 
iors are marked by arrows and the continuous lines are 
fits to the data of Eq. (fl"7|) . When CB is applied (inset), a 
cross-over between a regime with pj. « 1 to a decreasing 
power law is found. 

As a second characteristic of diluted networks, the 
emergence of correlations in the diluted network is dis- 
cussed. The theoretical result of Eq. (|3"3"j) suggests dis- 
assortative mixing in the diluted network in case central 
bias is applied. To observe these correlations in the sim- 
ulations, the mean nearest-neighbor degree is calculated 
as a function of the node degree. The result of our sim- 
ulation is shown in Fig. [3l As expected, no correlations 
are present in the original network (top curve in red, 
/ = 1). However, with 10% of its edges removed, the 
mean nearest-neighbor degree in the diluted net clearly 
decreases as the degree of the node increases (bottom 
curve in black). This is a clear indication of disassorta- 
tive mixing after sequentially removing a certain fraction 
of edges and this confirms our theoretical prediction. 
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FIG. 3: Mean nearest-neighbor degree k nn as a function of 
the degree of a node k, both for the original network (red 
line) and a network where only a fraction / = 0.9 of the links 
is present (black line). Clearly correlations which give rise to 
disassortative mixing emerge in the depreciated network. We 
used an original network constructed with the uncorrelated 
configuration model with 7 = 2.6, m = 2 and N = 10 J which 
is diluted using sequential biased percolation with a = 0.3, 
such that 7 = 3.3. The result was obtained with two network 
realizations, on both of them the percolation process was ap- 
plied four times. To reduce the noise level, the mean-neighbor 
degree is averaged over eight successive values. 



B. Simultaneous Approach 

1. Comparison with Theory and Sequential Approach 

This subsection deals with the iterated simultaneous 
approach as introduced at the end of Sect. Ill Bl To ex- 
amine the percolation transition, we search for the prob- 
ability to belong to the largest cluster, Too, as a function 
of the fraction / of included edges. Results are given in 
Figs. HO and [U 

For random edge removal it is shown in Fig. 0] that 
the simultaneous approach coincides with the sequential 
approach; in that case, only one iteration is necessary 
to attain / = 1. When, on the other hand, a > 0, the 
simultaneous approach can only be used up to a certain 
value of f u , smaller than one (See Eq. (|2"4"|) ). However, 
the iterative procedure can be used until all links are 
included. In general, we can include over 98% of the links 
with a relatively small number of weight-recalculations or 
iterations. For instance, in case a network with 7 = 2.5 
is subjected to percolation with a = 0.2, one finds that 
/„ = 0.51 and 15 iterative steps are necessary to reach 
/ = 0.98. 

The results of such iterative simultaneous approach can 
be found in figure [5J It is immediately clear that the se- 
quential and the simultaneous no longer coincide. This 
is not surprising since the definitions of edge retaining 
probabilities pij are different for the sequential and the 
simultaneous approaches when a 7^ 0. Furthermore, the 
sequential approach introduces correlations which are ab- 
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FIG. 4: The probability that a node belongs to the giant 
cluster, Voo, as a function of the fraction of retained links / 
for unbiased percolation (a = 0). We compare the sequential 
approach simulations (blue dots), the iterative simultaneous 
approach simulations (red dots) and the theory of generating 
functions (black line). We used an original network with 7 = 
2.5 and TV = 10 5 . The analytical results up to the fraction 
f u = 1 are obtained with the generating functions theory. All 
results are in very good mutual agreement. 
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FIG. 5: The probability that a node belongs to the giant 
cluster, Voo, as a function of the fraction of retained links 
/ for biased percolation with a = 0.25. We compare the 
sequential (blue line) and simultaneous (red dots) approach 
simulations. We used an original network with 7 = 2.5 and 
TV = 10 5 . 



sent in the simultaneous approach. Note, however, that 
the difference between the two approaches is not a mere 
consequence of the appearance of correlations in the se- 
quential approach. Indeed, if only the (dissassortative) 
correlations were present, the sequential approach should 
have a larger critical fraction f c than the simultaneous 
approach [6j; yet, we find the inverse to be true as ev- 
idenced in Fig. [SI At the point at which the weights 
are recalculated, the giant cluster probability Too of the 
iterative simultaneous approach shows kinks as a func- 
tion of /. For instance, a conspicuous kink appears at 
/„ = 0.69 in Fig. [5J The iterative simultaneous and 
sequential methods approach each other as / — > 1 and 
coincide at / = 1. 
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FIG. 6: The probability that a node belongs to the giant 
cluster, Voo , as a function of the fraction of retained links / for 
biased percolation with a = 0.25. We compare the sequential 
approach simulations (blue line), the iterative simultaneous 
approach simulations (red dots) and the theory of generating 
functions (black line). We used an original network with 7 = 
2.5 and N = 10 5 . The analytical results up to the fraction 
f u — 0.44 are obtained with the generating functions theory. 
This figure presents more detail of the critical region of Fig. [S] 



Note that the critical fraction f c which can be ex- 
tracted from Figs. [4] [5] and [6] is nonzero although it is 
expected to be zero for 2 < 7 < 3. This is a consequence 
of the finite-size effects which cause f c to scale with the 
system size according to Eq. (|3T))) (see also Fig. [1]). 

Although there are differences between the sequen- 
tial and the simultaneous approach, the differences are 
clearly not very large. The largest relative deviations oc- 
cur around f c while the largest absolute deviations are 
situated around the lowest /„ and are typically not more 
than 10%. Although we find the critical fraction f c of 
the simultaneous approach to be always larger than the 
one of the sequential approach, both values deviate by 
less than 10%. We conclude that the simultaneous and 
the sequential approach do differ, but the differences are 
not large and both approaches are qualitatively similar. 

Analytical results for the probability of the largest clus- 
ter Vqo in the regime f < f u can be obtained by solving 
Eq. (|45ap numerically for -ffi(l) which is then introduced 
in Eq. (|45b|) . For random edge removal (a = 0), /« = 1 
and thus a theoretical result is available for all /-values. 
Moreover, also for other values of a the generating func- 
tions theory calculation (see black line in Fig. [5]) agrees 
well with simulation results throughout the entire recon- 
struction process. The theoretical model is thus justified 
by the simulations. 



2. Properties of the Diluted Network 

We end the section with an overview of the properties 
of the diluted network. The focus lies again on the ap- 
pearance of correlations and on the cross-over. In the 
simultaneous approach (with / < /„), no cross-over can 




k 

FIG. 7: The probability pk to retain a node after depreciation 
(see Eq. (JT5J ) as a function of its original degree k for different 
values of / and using the iterated simultaneous approach. We 
used a network with 7 = 2.5, a — 0.5 and N = 10 and 
averaged over ten network realizations. Results are averages 
over ten percolation simulations. Since f u — 0.44, the first 
regime can be reached in a single sweep and thus no cross- 
over emerges for / = 0.1. The cross-over (indicated by arrow) 
becomes apparent for / = 0.9 when ten iterative steps are 
performed. This figure should be compared with the inset of 
Fig. [2] where the sequential approach was used. 



appear for the node retaining probability pk ■ Indeed, ex- 
pression (f2"2"| is exact and predicts a decreasing power law 
for pk, that is when / < f u . However, this expression is 
only valid as long as no recalculation of the weights is per- 
formed. As soon as we iterate the simultaneous approach, 
the results for the simultaneous and sequential methods 
start approaching each other. Since the sequential ap- 
proach contains a cross-over, we expect the appearance 
of a cross-over in the iterative simultaneous approach. 
Analoguous arguments apply to the appearance of corre- 
lations in the diluted network. 

The emergence of a cross-over in the iterative simul- 
taneous approach is indeed found and shown in Fig. [7] 
When 10% of the edges are included, no cross-over at all 
appears. This is consistent with our arguments since we 
can simply include this fraction of edges in one sweep. 
Also in the second sweep, no cross-over appears. How- 
ever, after 10 recalculations of the weights, the cross-over 
is undoubtedly present. Once again, our theoretical ar- 
guments are verified. Furthermore, the appearance of 
correlations is confirmed by our simulation results as ev- 
idenced in Fig. [5] Indeed, disassortative correlations are 
apparent only for large values of /, after several iterations 
have been performed. Note also that these correlations 
disappear only very slowly upon approach of the point 
/ = 1 where no correlations are present (see Fig. [3]) . 
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FIG. 8: Emergence of correlations in the diluted network in 
the iterative simultaneous approach. This figure shows the 
mean nearest-neighbor degree in the diluted network. We 
used an original network constructed with the uncorrelated 
configuration model with 7 = 2.6, m — 2 and N — 10 5 which 
is diluted using the iterative simultaneous approach to biased 
percolation with a — 0.3045, such that 7 = 3.3. The diluted 
networks contain, from bottom to top, 40 (red), 50 (black), 
90 (green), 95 (blue), 98 (purple) and 99.5 (grey) percent of 
the edges of the original network. To obtain such a network, 
resp. 1, 2, 5, 8 and 12 iterations have been made. The result 
was obtained with 10 network realizations, on all of them the 
percolation process was applied 10 times. To reduce the noise 
level, the mean nearest-neighbor degree was averaged over 8 
successive values. 



VIII. CONCLUSIONS 

We have performed a detailed study of biased percola- 
tion on scale-free networks with degree exponent 7 (with 
7 > 2) and shown that it is possible to tune a robust 
network fragile and vice versa. Biased percolation in- 
volves degree-dependent removal of edges, more specifi- 
cally, we assumed that the probability to retain an edge 
is proportional to (kikj)~ a with fej and kj the degrees of 
the attached nodes. For a > the bias is central since 
links between highly connected nodes are preferentially 
depreciated, while the converse, peripheral bias, corre- 
sponds to a < 0. Our most important result is that, at 
percolation, the properties of a network with bias expo- 
nent a and degree exponent 7 are the same as those of 
a network with bias exponent zero and degree exponent 
7 = (7— a) /(l — a), or degree exponent 7, depending on 
the sign and the range of a. Let us first elaborate on this 
main result, in the light of the present work and reca- 
pitulating arguments presented in the preliminary report 

For a > (with the restriction a < 1) the new de- 
gree exponent 7 > 7 governs the critical properties of 
the network that results when the percolation thresh- 
old is reached after biased depreciation. The exponent 
7 controls the large-degree behavior of the new degree 
distribution. This new P(k) is not simply scale-free but 



asymptotically scale-free. There is a cross-over value k x , 
so that for k < k x the exponent 7 is dominant and for 
k > k x the exponent 7 takes over. For 7 > 3 the biased 
depreciation process will reach the percolation threshold 
at a finite fraction of retained edges. The network is 
then fragile under central bias, regardless of whether the 
network is fragile (7 > 3) or robust (7 < 3) under ran- 
dom removal. For 7 < 3 the biased depreciation will (in 
an infinite system) not reach a percolation point since 
the critical fraction of retained edges, f c , is zero. The 
network, which is robust for random removal, remains 
robust under centrally biased removal. However, for a 
finite system f c is small but finite and scales with system 
size in a manner governed by the exponent 7, whereas the 
scaling properties of f c for random removal are governed 
by 7 < 7. We have shown, by analytic proof, that 7 gov- 
erns the percolation critical behavior for the case a > 0. 
Also our numerical results support this conclusion. 

For a < (with the restriction 2 — 7 < a), the criti- 
cal behavior at percolation is more subtle. Peripherally 
biased removal is less destructive than random deprecia- 
tion and it is possible that a network that is fragile un- 
der random failure becomes robust when peripheral bias 
is applied. Noting that 7 < 7 it is obvious that robust- 
ness is preserved for networks with 7 < 3. Conversely, 
fragility persists for sure when 3 < 7. However, it is not 
obvious what to expect when 7 < 3 < 7. The behavior of 
the new degree distribution P(k) for k > k x is, for a < 0, 
controlled by the exponent 7, so it would seem that the 
properties of the network under random failure are sim- 
ply not affected by peripheral bias. However, a finite-size 
scaling analysis at criticality reveals that the cross-over 
value k x is larger than the maximal degree in the net- 
work, implying that the new degree distribution will be 
controlled by 7 instead of 7, provided 2 — 7<a<3 — 7 
(we assume 7 > 3 since this discussion only makes sense 
for networks fragile under random failure). We conclude 
that sufficiently strong peripheral bias can turn a frag- 
ile network robust, and numerical evidence supports this 
conclusion. On the other hand, for 3—7 < a < 0, it is not 
clear whether the network stays fragile under peripher- 
ally biased failure when 7 drops below 3, which happens 
for a < (3 — 7)/2. This problem is still largely open to 
future investigation. 

Two distinct approaches by means of which a network 
can be reconstructed in a degree-dependent manner, the 
sequential and the simultaneous approach, have been in- 
troduced to perform the edge removal process. For the 
sequential approach, we obtained a very useful analytic 
approximation to the marginal distribution pk, which is 
the mean probability that an edge connected to a node 
with degree k is present in the network after reconstruc- 
tion. This analytic form clearly features the cross-over 
value k x which plays a crucial role in the network prop- 
erties. The simultaneous approach, which is a simpler 
scheme useful for a > and for edge number fractions 
below a value dependent on 7 and a, can be iterated so as 
to provide an alternative to the sequential approach (for 
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a > 0). The iterations introduce a history-dependence 
and lead to the emergence of fc x , rendering both recon- 
struction methods qualitatively similar. 

For both approaches the new degree distributions 
have been calculated and the degree-degree correlations 
emerging in the depreciated network have been charac- 
terized, by means of standard combinatorial methods. 
The main finding as regards the correlations is that the 
sequential approach causes disassortative mixing in the 
depreciated network when a > 0. 

For the simultaneous reconstruction approach, the ex- 
act (since correlation-free) percolation threshold f c is de- 
rived for central bias (a > 0) as a function of (non- 
integer) moments of the degree distribution, for 7 > 3. 
On the other hand, for 7 < 3 the exact finite-size scal- 
ing law for the vanishing of f c is obtained. These results 
fully demonstrate the validity of our exponent mapping 
© for central bias. 

A generating functions approach is introduced for 
degree-dependent edge percolation, extending previous 
work on random percolation. This approach allows to 
obtain the size distribution of finite clusters close to the 
percolation transition as well as other critical properties. 
If the network can be treated as a tree, which is valid for 
all finite clusters, the generating functions satisfy self- 
consistency equations. We have derived the extensions 
of these equations for degree-dependent percolation, al- 
lowing for correlated networks, and have shown that they 
reduce to the original equations provided no correlations 
are present. We have also derived the criterion for the 
percolation threshold for degree-dependent percolation 
and have shown that it reduces to the familiar Molloy- 
Reed criterion when correlations are absent. Further, 
our self-consistency equations reduce, for random perco- 
lation, to equations frequently encountered in the litera- 
ture. Our generating functions formalism is new in the 
sense that it extends known results on random percola- 
tion to biased percolation which may involve correlated 
networks. In the following, however, we draw further 
conclusions for the statistical properties (including criti- 
cal exponents) of uncorrelated networks only. 

Using the equivalence between the q — > 1 limit of the 
Potts model and edge percolation, we have shown that 
critical exponents for our biased percolation problem can 
be obtained from the Potts model free energy by extend- 
ing this equivalence to inhomogeneous (edge-dependent) 
couplings in the Potts model and edge-dependent removal 
probabilities in percolation. The generating functions 
approach has been combined with the extension of the 
Fortuin-Kasteleyn construction for the Potts model and 
with finite-size scaling in order to extract the critical ex- 
ponents of the percolation transition, for uncorrelated 
networks. We have found that the critical exponents 
are functions of 7, assuming that the degree distribu- 
tion after depreciation is governed by degree exponent 7, 
asymptotically for large degree. For 7 we obtain critical 
exponents that reduce to literature values of random per- 
colation simply by substituting 7 — > 7. However, in the 



more delicate regime 2 < 7 < 3 this correspondence is 
not satisfied. A critical assessment of this discrepancy is 
not given here, but left to future scrutiny. We conclude 
that, in all cases, the only way in which the bias expo- 
nent a enters in the critical exponents of the percolation 
transition, is through the new degree exponent 7. 

Furthermore, we have used numerical simulations to 
study the properties of the network after depreciation 
and near the percolation transition. We verified that 
robust networks can turn fragile under centrally biased 
failure and that fragile networks can turn robust under 
(sufficiently) strong peripherally biased failure, using the 
sequential approach. Although correlations are intro- 
duced in this approach, the results agree well with the 
predictions for uncorrelated networks. Also the cross- 
over behavior of the new degree distributions was tested 
and found to agree well with the analytical expectations. 
As regards correlations introduced by the sequential ap- 
proach, we have been able to verify the occurrence of 
disassortative mixing predicted theoretically for a > 0. 

The critical properties at percolation were checked by 
simulations using the (iterated) simultaneous approach 
and also compared with results obtained by simulations 
using the sequential approach. Specifically, we have 
found that for biased percolation the sequential and the 
(iterated) simultaneous approach give rise to different re- 
sults. In particular, the size of the giant cluster pre- 
dicted by the generating functions theory agrees very 
well with the simulations for the (iterated) simultaneous 
approach. Nevertheless, the differences are often small 
and we may conclude that both methods are qualita- 
tively similar. Finally, we have also provided evidence 
for the theoretically expected appearance of cross-over 
effects and degree-degree correlations for the (iterated) 
simultaneous approach. Overall, we conclude that good 
agreement has been found between simulations and the- 
ory. 



IX. APPENDIX 

There exists a more formal way to prove that the de- 
gree distribution in the diluted net satisfies P(k) cx k ~ 7 
for large degrees k when central bias is applied. 

We start from the degree distribution in the diluted 
network P(k) which was calculated in Eq. (|26|) . For large 
values of k, i.e., k 3> k x , and a > 0, the probability of 
retaining a node of degree k falls off as pk oc k~ a : this 
is valid using the sequential approach (see Eq. (|2"D|0 . as 
well as the simultaneous one (see Eq. (|2"2l ). 

If both kpk and fc(l — pk) are large, the binomial dis- 
tribution can be approximated by a normal distribution 
with mean kp^ and variance kpk(l — Pk)- The latter con- 
dition is always true if we apply CB, since then l-pt«l 
for large k as edges between the most connected nodes 
are almost certainly removed. Since kpu oc fc 1- ", the first 
requirement holds only if a < 1. 

Inserting the normal distribution with variance kp^il— 
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p k ) w Cok 1 a , with Co a constant, in Eq. (|26|) and ap- 
proximating the sum by an integral yields 

P( fc) oc^ rffcfc -^exp^- l 2Cofcl _J J- 

(80) 

Now we introduce the auxiliary variable u = k/k ^ ^ 
and rewrite the integral as follows: 

l-T 1 f°° „_! 

P(fc)ocfc TT ^ T 2 /__ B 
A 1-° 

XeXP (- fc [ 2C ^-° J" (81) 

For k — > oo, the integrand has only non- vanishing values 
in a neighborhood Aw « 1/Vk around u c = C ~ X '^ a '. 
If < a < 1, the lower bound of integration vanishes 
for large fc. Thus u c certainly lies in the domain of in- 
tegration and the integral can simply be approximated 
by Ci/Vk with C\ a constant. After some trivial power 



counting, we arrive at 

P(fc)ocfc-^. (82) 

Thus we obtain the anticipated behavior for CB. The 
exponent 7 controls the decay of the degree distribution 
in the diluted network at large k. 

The case a = 1 can be studied using a Poisson dis- 
tribution approximation for the binomial factors, which 
results in P(k) being a Poisson-type degree distribution 

P{k) ~ k ~ jk\ for large k. Thus a < 1 is a natural 
restriction, because the scale-free behavior is destroyed if 
stronger CB is applied. 

We still have to examine the situation for a = 0. 
Then the lower bound of the integral becomes 1, while 
u c = I /Co- The constant Co is nothing but the fraction 
of edges preserved after the depreciation process. Thus 
Co < 1 and u c also lies in the integration domain. Thus 
the previous arguments apply as well to the random node 

removal process. We conclude that indeed P(k) ~ k 7 
for < a < 1 which proves the intuitive conjecture given 
in the Introduction. 
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